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Resonance diffraction in the periodic array of graphene micro-ribbons is theoretically studied 
following a recent experiment [L. Ju et al, Nature Nanotech. 6, 630 (2011)]. Systematic studies over 
a wide range of parameters are presented. It is shown that a much richer resonant picture would be 
observable for higher relaxation times of charge carriers: more resonances appear and transmission 
can be totally suppressed. The comparison with the absorption cross-section of a single ribbon 
shows that the resonant features of the periodic array are associated with leaky plasmonic modes. 
The longest-wavelength resonance provides the highest visibility of the transmission dip and has the 
strongest spectral shift and broadening with respect to the single-ribbon resonance, due to collective 
effects. 

PACS numbers: 42.25.Bs, 41.20.Jb, 42.79.Ag, 78.66.Bz 



The ability of graphene to support electromag- 
netic waves coupled to charge carriers [graphene sur- 
face plasmons (GSPs)] is very interesting from the 
point of view of many physical phenomena related 
to surface plasmons (SPs) 12 . An additional inter- 
est is related to graphene's flexibility, sensitivity to 
external exposure and two-dimensionality (2D) that 
have a variety of possible applications GSPs have 
been intensively studied theoretically^^ i n gra phene 
sheets, and also in g ra p hene ribbon jH H 16 l 18 r p-n 
junctions^ and edges ^Hl and recently have been ob- 
served experimentallyP^EIl 

In metal films, the excitation of the SP modes had 
been experimentally and theoretically studied for pe- 
riodic ultrathin structures (< 10 nm-thick), both for 
arrays of slabs^HUl an( l arrays of holes and disks 26 28 . 
It has been shown that these systems present trans- 
mission peaks with high visibility (including total sup- 
pression of reflection) and absorption resonances. The 
natural continuation of this research was to check if 
this property could still hold for the 2D limit i.e., 
for a layer of one- atom thickness. Recently, experi- 
ments have shown that GSP resonances in a periodic 
array of graphene ribbons (PAGR) have remarkably 
large oscillator strengths, resulting in prominent room- 
temperature optical absorption peaks P4 

In this paper we present a theoretical study of the 
electromagnetic response of PAGRs, including absorp- 
tion, transmission and reflection coefficients. We con- 
sider both the parameters corresponding to the exper- 
iment and their variation over a wide range. Specifi- 
cally, we focus on the dependencies upon the relaxation 
times of charge carriers r and the width-to-period ratio 
(which in the experiment was fixed to be 1/2). We look 




FIG. 1: (Color online) The geometry of the studied system: 
a periodic array of graphene micro-ribbons of width W and 
period L, with a normally-incident electromagnetic wave 
having the magnetic field along the ribbons. The array 
is placed between two dielectric half-spaces with dielectric 
constants Si and €2- 

for the configurations in which GSP-induced absorp- 
tion is enhanced and where other GSP-assisted effects 
are much more pronounced. Our analysis can thus be 
used for further efficient observation of GSPs and their 
use for applications e.g., ultra-thin voltage-controllable 
THz absorbers. 

Figure [I] schematically represents the periodic ar- 
ray of graphene ribbons under study. The PAGR is 
located at z = and is illuminated by a normal- 
incident monochromatic plane wave (having vacuum 
wavelength A), with electric field pointing along the 
x— direction. The period of the PAGR is L, the width 
of the ribbon is W and the dielectric permittivities 
of the superstrate and substrate are e± and £2, re- 
spectively. The graphene ribbons are modeled using a 
2D conductivit y a, c omputed within the random-phase 
approximationPSE] Room temperature, T = 300K is 
considered throughout the paper. 

Due to diffraction, the PAGR generates an infinite 
discrete set of plane waves n G Z with ^-components of 
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FIG. 2: (Color online) Absorption and gate-induced 
change of transmission spectra for PAGR with the di- 
electric cladding and relaxation time corresponding to the 
experiment!^ e\ — 3, £2 — 4, r = 0.25 ps. (a) shows the 
absorption spectra A for different values of chemical poten- 
tial \i. The continuous curves are for the period L — 8/xm, 
while the discontinuous ones are for L — 4/xm. For both 
cases W — L/2. In the inset the relative change of trans- 
mission with respect to the sample at charge neutral point 
ST = — (T — Tcnp) /Tcnp is shown for the same L and 
W as in the main figure. In (b) the spectra for A (contin- 
uous curves) and ST (discontinuous curves) are shown for 
different widths of the ribbons W for L = 8/xm, \i — 0.2eV. 
The case of W = 4/xm is also shown for the free-standing 
graphene [si — £2 — !)• 



the wavevectors k nx = nG, G = 2tt/L being the short- 
est vector of the reciprocal lattice. The fields in the 
dielectric half-spaces can then been presented in the 
standard form of Fourier-Floquet expansion. Match- 
ing the fields at the interface z = results in an infinite 
set of linear equations for the amplitudes of diffracted 
waves. The direct calculation of the diffraction ampli- 
tudes in the truncated linear system is simple to imple- 
ment and, additionally, provides qualitative informa- 
tion on spectra of GSPs. However, the convergency of 
this procedure with respect to higher harmonic consid- 
ered, TV, is poor for the chosen polarization. Therefore, 
for each geometry and for the lowest scattering time 



r considered, the modal expansion calculations have 
been checked by the finite elements method (FEM) 
realized in COMSOL. Once the value of N needed to 
achieve convergency is found, the faster modal expan- 
sion method can be used to study the dependency with 
r of the scattering coefficients. 

Let us start our analysis of the electromagnetic prop- 
erties of PAGR by the geometry considered in the ex- 
periment described in Ref.|20| For this, we take e\ — 3 
(ion gel) and £2=4 (Si02), ignoring thus the effect of 
a finite thickness of the Si02 layer and possible related 
Fabry- Perot type resonances. The scattering rate is as- 
sumed to be 4THz (r = 0.25ps). As in the experiment, 
the transmission coefficient is compared with the one 
at the "charge neutral point" Tcnp, where the chem- 
ical potential is very small (we take /i = 10 _2 eV). 

First of all, following the experimental study, we 
consider the variation of the spectra with the change of 
the chemical potential \i and the period L, for a fixed 
ratio, W/L = 1/2. Fig. [2] (a) shows the transmission 
coefficient change ST = — (T — Tcnp) /Tcnp an d ab- 
sorption coefficient A as a function of the wavelength. 
For each value of L and \i there is a resonant maximum 
in both A and ST spectra. As will be seen below, the 
resonance is related to the excitation of the longest- 
wavelength GSP in each ribbon. The maximum reso- 
nant absorption increases with the increase of doping, 
due to both the resonance shift to a less absorptive 
frequency region and to the higher number of charge 
carriers that get involved in the plasmonic oscillation. 
In accordance with the experiment, this resonance 
blue-shifts when either \i increases or L decreases. 
This behavior can be explained from the condition 
for GSP resonance in the ribbon, which approximately 
satisfies 16 W ~ nXcsp/Z- Here Xgsp — 27r/Re(/cGS , p) 
is the GSP wavelength, and n measures the number of 
half- wavelengths that fit within the ribbon width for a 
certain mode. In the considered frequency range, the 
intra-band Drude-like term dominates in the conduc- 
tance, so Re(kcsp) — hu 2 /(2ao[ic), where ao is the 
fine-structure constant. Substituting Re(kcsp) into 
the resonance condition, we have for the resonance 
wavelength \ res ~ ^2nchW/ (nao/n) ex y/W/fi. 

Further insight into the absorption process can be 
gained from going beyond the W/L = 1/2 ratio con- 
sidered in the experimentPSl. Figure J2^b) shows the 
spectra for A and ST for different widths of the ribbons 
W, at the fixed period L = 8/im and for \i = 0.2eV. 
For larger values of W the resonance shifts to longer 
wavelengths, where graphene is more absorptive and, 
correspondingly, the peak broadens. Interestingly, the 
propagation length of GSPs increases at longer wave- 
lengths since the increase of Re(cr) is overcompensated 
by the lower GSP confinement. Nevertheless and de- 
spite the presence of resonances, for the considered 
r = 0.25ps, the maximum of A grows with W reaching 
its maximum for a continuous graphene sheet W = L 



3 



10" 



10" 



10 





- 1 t 

3 2 


t 

1 




W = 0.\L 




— • 40 ps single ribbon 
40 ps sheet 
10 ps sheet 
0.25 ps sheet 











0.25 ps 






10 ps 






40 ps 










l\ R 




■ 4 \ 





X, \im 



10 



X, \im 



10 



FIG. 3: (Color online) Absorption A [in panels (a),(c),(e)], reflection R and transmission T [in panels (b),(d),(f)] spectra 
for free-standing periodic array of graphene ribbons with different values of ribbon width W and relaxation times r. Panels 
(a)-(b), (c)-(d), (e)-(f) correspond to W = 0.8, 4 and 7.2/xm, respectively; in all cases the period is L = 8/xm. The inset 
to (b) represents the colorplots for the electric field modulus |E| and the real part of the ^/-component of the electric field 
E y in the vicinity of the ribbon. The colorplots are marked with respect to the numeration of the peaks in panel (a). 
The inset to (d) shows the dependencies of the maximum values of A and R and the minimum value of T as a function 
of r, for longest- wavelength resonance shown in (c)-(d). The horizontal discontinuous line in the inset sets the maximal 
possible value of A. The dotted lines in panels (a),(c),(e) represent the absorption by a continuous graphene sheet. 



(this dependence is almost linear, as shown by cal- 
culations for intermediate values of W not presented 
here). In other words, for small values of r the ef- 
fect of GSP-induced absorption is weak, so that the 
absorption is approximately proportional to the area 
covered by graphene. Actually, as will be rendered 
in Fig [3^c) , the computed absorption is higher for a 
continuous graphene sheet than for a PAGR with the 
W/L = 1/2 ratio considered in the experiment. 

Our calculations show that the GSP-absorption 
effect would be greatly enhanced for higher val- 
ues of r, which are currently associated to free- 
standing graph ene s heets and their much higher elec- 
tron mobilities ! 33 * 34 * In order to differentiate between 
the effects of changing the dielectric environment and 
changing the relaxation time, Fig. j2^b) renders A and 
5T for a free standing PAGR, with W/L = 1/2 and 
r = 0.25ps. The resonance in the free-standing PAGR 
blue-shifts and is get narrower than that of the corre- 
sponding PAGR with dielectric surrounding (which is 



related to the shift of GSP dispersion curves), but all 
the discussed tendencies with the change of \i and W 
are the same. Similarly and even though the maximum 
absorption in the PAGR has increased, the absorption 
in the spectral window considered is below the one for 
a continuous graphene sheet. 

This situation changes at higher relaxation times. 
Figure [3] illustrates the absorption, transmission and 
reflection spectra for the suspended PAGR of different 
relaxation times r = 0.25, 10 and 40ps and different 
widths of the ribbons W = 0.1,0.5 and 0.9L, for the 
period L = 8/im. To make a set of the resonance 
peaks more visible (specially those appearing at lower 
wavelengths) the absorption is presented in logarithmic 
scale. For each ribbon width the absorption coefficient 
is compared with the absorption cross-section (Asr) 
for a single ribbon of the same width at r = 40ps. For 
better comparison, the Asr is normalized so that its 
value coincides with A at the shortest wavelength in 
the considered spectral interval. 
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Each peak on the absorption spectra corresponds to 
a GSP resonance in the ribbon. Increasing W increases 
the number of resonances that appear in the spectral 
window considered. These resonances correspond to 
the excitation of either GSP waveguide- or edge- type 
modes with zero value of ^-vector in the ^/-direction. 
These are leaky modes, resulting from the GSPs dis- 
cussed in Ref. 16, with the prolongation of the dis- 
persion curves inside the light cone down to the value 
Re(k y ) = 0. The field distribution around a ribbon 
corresponding to the last three peaks in the absorp- 
tion spectra is shown in the inset of Fig. [3^a). The 
two highest-wavelength modes result from the degen- 
erate edge GSPs while the rest of the resonances corre- 
spond to excited waveguide-type GSPsP^ Notice that, 
in practically all cases, the absorption spectra for the 
array and the single ribbon are approximately equal 
(independently on the value of r) . The only exception 
occurs for the resonance appearing at the longest wave- 
length, and only for narrow gaps between the ribbons 
(gap width < 0.2L), when the GSPs in neighboring 
ribbons hybridize. 

In the symmetric dielectric environment considered, 
A in the graphene array can not exceed its maximal 
value^ Am = 1/2. But importantly, even for small 
ratios W/L, for sufficiently large values of r the GSP- 
induced absorption in PAGR can not only be higher 
than the absorption corresponding to lower r, but can 
also largely exceed the absorption in the continuous 
graphene sheet [see Fig. |3ja),(c),(e)]. 

The calculations rendered in Fig. [3jb),(d),(f) show 
that absorption peaks are complemented by peaks in 
reflection R and dips in transmission T, with the 
longest-wavelength resonance presenting the deepest 
minimum in T. 

Let us now focus on the longest-wavelength reso- 



nance at W/L = 1/2 [see Fig.^c) and Fig.^d)]. The 
inset to Fig.|3jd) renders the dependency on relaxation 
time of the maximum values of A and R together with 
the minimal value of T. As seen, the evolution of A max 
with r is not monotonous, with A max reaching the 
optimum value at r c± 3ps. Conversely, T min mono- 
tonically decreases with r, while R max monotonically 
increases with it. Importantly, the minimal value of 
the transmission in the resonance keeps its low value 
T min < 10% down to r > 7ps, and T min < 2% for 
r > 30ps. Taking into ac count hight values of mobili- 
ties in suspended samples j 33 l 34 l these deep transmission 
minima could be observed experimentally. 

To conclude, we have studied the transmission, re- 
flection and absorption resonance THz spectra in peri- 
odic arrays of graphene ribbons. The resonance effects 
are related to the leaky plasmonic modes existing in 
individual ribbons and the modes corresponding to dif- 
ferent ribbons are very weakly coupled to each other. 
The highest-wavelength resonance provides the max- 
imal visibility of the transmission dip and reflection 
peak, with its resonant character surviving even for 
the low relaxation times present in graphene samples 
on a substrate. As this mode is the less confined, it 
is the most strongly perturbed by the periodicity of 
the array. The samples with higher relaxation times 
allow for more resonances being visible and provide 
very deep transmission minima. We have shown that, 
in ribbon arrays with sufficiently high relaxation time, 
the absorption can be substantially higher than the 
absorption in the continuous graphene sheet. 
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